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Abstract 

We develop and apply an effective analytic theory of a non-collinear, broadband type-I para- 
metric down-conversion to study a coupling efficiency of the generated photon pairs into single 
mode optical fibers. We derive conditions necessary for highly efficient coupling for single and 
double type-I crystal producing polarization entangled states of light. We compare the obtained 
approximate analytic expressions with the exact numerical solutions and discuss the results for a 
case of BBO crystals. 
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I. INTRODUCTION 



Sources of polarization entangled photon pairs are important in a wide range of modern 
applications. Efficient generation of photon pairs led to an opportunity of testing sev- 
eral fundamental problems of quantum mechanics and to a practical realization of various 
quantum-information or cryptographic schemes. Among all accessible entanglement sources, 
down-conversion crystals turned out to be the most effective and easy to implement so far. 
Therefore the development of down-conversion techniques started to be an important issue 

n 

for many branches of quantum optics 1]. In a search for enhancement of counting rates, it 
has been noticed that focusing of the pump beam can essentially increase the efficiency of 

n 

the process |2| . Another important observation was that coupling of the generated photons 
into single mode fibers can optimize the collection efficiency and is a very practical method 
of preparing photons in well defined spatial modes 0|. This idea has already inspired some 
authors to study various aspects of fiber coupling of the down-converted photon pairs 
Efficient fiber coupling of the entangled photon pairs can also help in enhancing the fiber 
communication efficiency by encoding additional information in the polarization degree of 
freedom Q|. 

The paper delivers a detailed analysis of the fiber coupling efficiency and properties 
of the coupled type-1 down-conversion photon pairs in a general case of a pumping with 
ultrashort pulses. To our best knowledge this issue has not been analyzed so far. One of the 
advantages of using the pulsed pump sources is the fact, that generated photon pairs arrive 
at the detectors at a well-defined temporal window which may be very helpful in a number 
of experimental schemes ^. 

We develop an approximate analytic theory of type-1 down-conversion to study the de- 
pendence of the coincidence spectrum and coupling efficiency on several relevant parameters 
characterizing the system such as properties of the nonlinear medium, pumping beam, spec- 
tral filters, geometry of the setup, etc. We point out what are the optimum (for the coupling 
efficiency and a spectrum separability, which is necessary for a high-purity polarization en- 
tanglement generation [3]) settings of the setup and compare the result with the numerical 
calculations. 

The paper is organized as follows: in Sec. EJwe derive explicitly a wave function of the 
type-I down-conversion generated photon pairs and discuss the applied approximations. In 
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Sec. mil we discuss the coupling process, properties of the spectra of the coupled photons, 
derive the spectrum separability condition and derive expressions for the coupling prob- 
ability. In Sec. HVl we apply the obtained results to discuss the case of a double type-I 
crystal generating polarization entangled states and finally Sec. |V] concludes the paper and 
in the Appendix we derive a relation between geometrical parameters of the crystal and a 
dispersion relation. 



II. WAVE FUNCTION 



Consider a type-I down-conversion process taking place in an anisotropic medium with a 
second-order nonlinearity. Interaction between an extraordinarily polarized laser pump beam 
Ep{r,t) (which will be assumed to have a finite time) and the nonlinear crystal generates 
pairs of daughter photons whose properties are fully determined by the properties of the 
pump beam, indices of refraction and geometry of the process - see Figure ^ In the first 
order perturbation, only a single photon pair can be created in the down-conversion and the 
output quantum state has the approximate form: 

dt J^dhEAr,t)(^E^-\r,t)) |0), (1) 

where an electric nonlinear polarizability x is a parameter determining the interaction 
strength, V is volume of the crystal and a negative part of the electric field operator is 

, uj) (a is an annihilation operator of an ordinar- 
ily polarized mode, at ^ = 0). The integral extends to all positive frequencies u and all 
perpendicular wave- vectors k^, ky. Such a choice of the plane wave decomposition will turn 
out to be very handy, since the perpendicular wave-vector components and frequency do 
not change on dielectric media boundaries P, 0]. Let us introduce signal and idler mode 
vectors s = (fcs,a;, ks,y,ujs) and i = {h^x, k\^y,uji) to decompose the square of the electric field 
operator appearing in the expression (PJ). Then the 2;th components of the wave vectors are 
given by a dispersion relations inside the anisotropic medium k^^z = k°{s) and k[^z = 
(index "o" refers to ordinary wave of the considered anisotropic medium - see Appendix) 
and the output quantum state can be written down using a quantity \l/(s,i) interpreted as 
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a probability amplitude for a photon pair to be created in modes s and i: 



1^) ~ |0)-^| J d'sd='i^(s,i)at(s)at(i)|0). (2) 

The first order approximation that we will study is justified when x <C 1. The properties 
of the generated photon pairs are described by the function \l/(s, i) and this is a typical 
situation in many down-conversion experiments. 

Suppose that a type-I down-conversion medium occupies the area —^<z<^ and it is 
pumped by a light beam propagating towards z direction. Then the probability amplitude 
\l/(s, i) at the end of the crystal has the following form: 



OO POO POO 



^(s,i) = / dt I dx I dy I dzEp{r,t)x 

(3) 



OO J —OO J — OO 



-i[(ks+ki)-r-(t<;s+i^i)t] „i(fcs,z+M,z)|- 



where the last exponent is a propagation factor of the amplitude from z = to z = ^. 

Let us express the pump beam electric field in terms of its Fourier transform defined 
at z = -\: Ep{r,t) = J dfcp,^. dfcp^^ dcjp e*(''p'''"'^p*)e*^P'^"i'Ep(fcp_a,, /cp^^, cjp), where we have 
introduced vector p = {kp^x, kp^y,^p) parameterizing the plane wave decomposition of the 
pump beam. To calculate the probability amplitude (IHl) we need the dispersion relation 
kp,z = kl{p), which is well known for uniaxial crystals - see Appendix. Inserting this 
into the equation and performing all the integrals yields 



0: 



^'(s,i) = Ep(s + i)Lsinc |^^5_(s,i)^ exp (^i^6+is,i) ] , (4) 
where 5± {5- is a so-called phase mismatch parameter) equal: 



5±(s,i)=fc:(s + i)±fc:(s)±fc:(i). (5) 

From the form of equation (jlj) it is clear that the process occurs most effectively in the 
directions for which the longitudinal wave vector components meet the phase matching 
condition 5_ = 0. To determine this condition explicitly one has to use Sellmeier expressions 
for the extraordinary and ordinary index of refraction as a function of the vectors s and i. 
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Since the expression ^ is not analytically integrable, which will turn out to be crucial 
in our analysis and the explicit form of the mismatch parameter 6- is rather complicated, 
we will perform certain approximations before starting any further analysis. First of all, 
let us perform a Gaussian approximation of the sine a; function appearing in with the 
expression e~^'^'^. This rough approximation is fair as long as we consider only photons 
generated near the phase-matching region of wave vectors and frequencies where the down- 
conversion occurs most effectively. In this narrow range of parameters, i.e. when <^ 1 
the wave tails of the sine x function give insignificant contribution to the overall probability 
distribution of generating a pair of photons. 

Secondly, we will expand 6± to the first order in the Taylor series around a point for 
which the phase-matching condition 5_ = is met in the horizontal plane. If we fix the x 
axis to be horizontal and the y axis to be vertical, our expansion takes place around the 
point So = (^^^,0,u}q) and io = (—^^,0,070), where cuq is the phase-matching degenerate 
frequency and 6*0 is the angle of propagation of the degenerate daughter photons relative to 
the z axis outside the crystal. These two parameters can be determined from the properties 
of the down-conversion medium (given the Sellmeier formulas for extraordinary indices of 
refraction and direction of the crystal axis, one can calculate c^o and 6q directly from the 
condition (5_(so,io) = and the definition (jS))). We obtain an approximate form of the 
probability amplitude (jH): 

^(s, i) ^ Ep{s + i)Lexp 

where we have introduced the following notation: Vf{xi,...XN) = JZiLi 'H'^^i- ap- 
parent that in our approximation all the relevant information about the structure of the 
nonlinear medium is contained in the derivatives of 6±. Particular components of these 
vectors depend on the properties of both, extraordinary and ordinary indices of refraction, 
as well as on the spatial orientation of the optical axis. 

For the reasons that will be made clear in the further parts of the analysis we shall choose 
the optical axis of the crystal in a plane oriented at the angle 45° in respect to the horizontal 
plane. In this case all the derivatives considered can be approximately expressed with only 
4 real parameters (moreover, it will turn out that one of them is irrelevant, so we end up 
with only 3 parameters) - see Appendix: 
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iV5. 



(6) 
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dS± d6± 7 

"7- 



d5± dS± 



(7) 

where 9'q is the angle of propagation of the degenerate daughter photons relative to the 
z axis inside the crystal proportional to the 9q (proportionality constant is given by the 
ordinary index of refraction), 7 (equal to the derivative of in the direction of projection 
of the optical axis onto the plane transversal respect to the z axis) is a walk-off angle of the 
extraordinary ray of frequency 2ujq incident in the z direction and A/^-t^^ = (31{2ujq) -^(3"{ijJq) 
are the combinations of the extraordinary and ordinary inverse group velocities towards 
z axis for the frequencies 2uoq and cuq, respectively. We have also assumed that the z-ih 
component of the ordinary inverse group velocity I3°{uj) weakly depends on the direction of 
propagation. 



III. COUPLING PHOTON PAIRS INTO FIBERS 



As we have seen, the full information about the modal structure of the generated photon 
pairs is contained in the function ^'(s, i) - the probability amplitude of emission of a photon 
pair of certain frequencies towards given two directions. In the next step of our analysis, 
we will study how these photon pairs are coupled into single-mode optical fibers. We may 
think of such fibers as a pair of spatial filters letting through the photons occupying only 
selected transverse modes u{x,y,u!) propagating through the fiber. 

Consider a scheme in which photons generated in the crystal pass through interference 
filters and lenses coupling incident light into the single-mode fibers. Let the fibers be oriented 
in the horizontal plane at the angles and 6i with respect to the z axis - see Figure H The 
coupled modes can be easily characterized if we consider an inverse problem: what is the 
structure of a light mode emerging from the fiber? The output mode u{x,y,u!) can be 
approximated by a Gaussian beam incident on the down-conversion crystal. Let us consider 
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a setup in which the coupled beam propagating in the horizontal plane (common with 
the pump beam) crosses the output face of the crystal of a distance h from the point of 
intersection of the pump beam and the crystal surface - see Figure ^ In the simplest case, 
when the output beam is focused on the surface of the crystal, it is parameterized only 
by the angle of incidence 9^-, the distance h from the pump beam and the waist w. The 
Fourier-transformed mode function at the output surface of the crystal for the signal arm 
(and analogously, only with —h instead of h for the idler arm) has the form: 



w 



exp 



w 



+ k. 



s,y 



UJs9s 



(8) 



where Ap^uj) = exp ^— ^'^^a'^^^ ) characterizes transmissivity of the interference filters charac- 
terized by their spectral bandwidth ap. The probability amplitude iIj{uJs,uJ]) for the photon 
pair of frequencies Ug and Ui to be coupled into the fibers reads: 



ijjius^uji) = J d/Cs,xdfcs,s/dfci,^dfci,y^(s, i)u*(s)M*(i). (9) 

Notice, that the above expression involves integrating of the probability amplitudes, not 
probabilities, because photons with different wave vectors, but the same frequency are in- 
distinguishable after coupling into the fibers. 

To write down explicitly the amplitude ipl^s, t^i) it is useful to introduce a set of effective 
parameters describing the process: a relative pump beam walk-off parameter F, a relative 
shift of the degenerate photons and an effective crystal length C: 



e 



LO'o 



w 



(10) 



v/1 + 02/2 + FV2 
These parameters appear naturally in the below expressions. 

To proceed with the analysis, let us assume that the pump field is a Gaussian beam 
focused on the crystal so that the Rayleigh range is much longer than the crystal length: 



Ep{kx, ky,uj) 



Ap{uj)wpexp {kl + ky) 



where the pump beam spectrum is Ap{uj) 
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exp ^— ''^^^2°'' j and the fibers are oriented at the angles for which the degenerate photons 
are collected the most efficiently: 9s = Oq and Oi = —9q. 

The explicit form of the probability amplitude modulus |^/'(a;s,a;i)| reads: 



|V'(^s,^i)| 



STlWpC 



exp 



1 



2 + r2 2 



exp 



2(2 + r^)/^;=^ / _ LO'^ 1 + y 
" 2 + e2 + r2 \ 2 i + r72y' 



(11) 



where uj — {cUg — uJojCUi — cuq) and particular components of the appearing spectrum matrix 



n 



^ss ^si 
^is ^ii 



read: 



(12) 



+ 



1 1 



C \/2c + 2^2 y 2c2(w2 + 2w2) ^2 



+ 



c V2c ^w^ + 2wl ) ' 2c2(w2 + 2w2) al 



^si — ^is 



OoP'o wT/e 
V2c ^^2 + 2w 




w'^WpOl 



+ 



2c2(w2 + 2w2) a- 



(13) 



Let us notice, that the diagonal elements are not exactly the same. The symmetry of the 
coincidence spectrum is broken by the non- vertical orientation of the crystal's optical axis. 
However, for small angles ^0 the spectrum matrix f2 simplifies and the asymmetry of the 
spectrum disappears: 



£2A/^ 



2„,,2/i2 



+ 



WW 



P"0 



2c2(w2 + 2w2' 



£2A/32 



w'^WpO^ 



2c2(w2 + 2w2) ^ Ct2' 



(14) 
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Non-separability of the function ilj{uJs,uJi) is a signature of frequency entanglement, a fea- 
ture that is undesirable if we are interested in polarization entanglement. To minimize this 
entanglement one has to adjust the parameters so that the off-diagonal elements of the spec- 
trum matrix Q become relatively small, i.e. when ^ reaches its minimum. Inspection of 
expressions (fT^ shows that for small angles and small mode waists w and wp the separa- 
bility can be accomplished only with the use of narrow-band interference filters. However, 
for large mode waists one can separate the spectrum by fulfilling the following condition: 



'4) 8 ' ~ al 



(15) 



2c2(w2 + 2wl) 8 (tI ■ 

Notice that this condition does not depend on the interference filters' bandwidth ap. 

For concreteness, we have analyzed the coincidence spectrum (fTT|) on the example of a 
BBO crystal cut for the degenerate frequency ujq = 780nm and the emission angle 6q = 1.4° 
(using the Sellmeier formulas for the indices of refraction [12] one can calculate that in this 
case A(3_^z = 2.06 ■ 10^^), pump beam spectrum width 5nm to and the interference filter's 
width 17nm. The results are shown on Figure El for several mode waist diameters. We see 
that when the beams are strongly focused, the spectrum \ip{ujs,uj[)\'^ becomes inseparable, 
and the separability can be achieved only for unfocused beams. In the considered case the 
optimum waists are approximately w = Wp ~ 2.5mm. 

Expression (fTTj) can be used to maximize the coincidence probability Q over all possible 
choices of the parameter h. The optimal shift h of the coupled beam in respect to the 
pumping beam at z = Z2 is: 

One can notice that in the absence of the walk-off effect affecting the propagation of the 
pumping beam (for F = 0), the optimal choice of h involves crossing of the pumping mode 
and two fiber coupled modes exactly in the middle of the nonlinear medium. The modifica- 
tion comes from the fact that the pumping beam is transversely shifted due to the walk-off 
effect during the propagation in the crystal. The effects of the walk-off become more signif- 
icant when all the modes are strongly focused on the crystal, because only the relative shift 
of the beam in comparison to its waist and the waist of the coupled modes is relevant. 
To obtain the total probability p of coupling any photon pair into the two fibers we need 
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to integrate the probability |'?/'(co's, over all possible frequencies. Notice that this time 
we integrate the probabilities, not their amplitudes as the photons of different frequencies 
are perfectly distinguishable: 



The above integration can be also done analytically if we extend the lower limits to — oo. 
In this case, for the optimum choice of h we obtain a relatively simple expression for the 
probability p: 



Before we continue with a detailed analysis of the above result, let us stress once more, that 
it is valid only when the Rayleigh range of all the considered beams is much longer than 



This condition guarantees that the diffraction of the pump beam and generated photon pairs 
inside the crystal can be neglected. 

Several conclusions about the efficiency of the coupling process can be drawn from the 
expression (fTHj) . First and the most obvious is that the overall probability p does not 
depend on A/J^ ^. Coupling efficiency is also a decreasing function of the angle 6*0 (and 
proportional as it appears only in the denominator of the expression (fTBj) . This can 
be easily understood, since increasing 6q (depending on the geometry of the crystal and its 
indices of refraction) decreases the spatial volume of an overlap between the pumping mode 
and coupled modes. Another clear conclusion is that the walk-off effect represented by F 
and responsible for the transverse shift of the pump beam during the propagation through 
the crystal, also diminishes the coupling probability. The next relevant parameter, inverse 
group velocity difference of the extraordinary and ordinary beams propagating towards z: 
Af3^^z suppresses the down-conversion process. This effect is also straightforward, as A/?_^^ 
is inversely proportional to the down-conversion spectrum width, so the larger AjS-^z is, the 
fewer frequencies are efficiently down-converted. 

The next issue of practical interest is how the probability p depends on the parameters 
adjustable for a given nonlinear material, namely the crystal length L and the mode waists 
w and wp. Again, we have analyzed this issue on the same example of the BBO crystal. 




(17) 




(18) 



the thickness of the crystal, i.e. for ^kQW^ ^ L, where kQ is a characteristic wave-vector. 
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For such setting we have calculated the pair coupling probability p given by Eq. (fTHj) and 
compared it with the result of a numerical calculation. In the figureElwe show the probability 
p as a function of the crystal length L for fixed w = wp = lOO/im (solid line) and an analogous 
numerical result (dashed line). The only approximation that the numerical calculus involves 
is the expansion of the 6± functions to the second order of the Taylor series. In the numerical 
case we keep the sine x function unchanged, as well as we do not assume that the Rayleigh 
range of the considered modes exceeds the crystal length. 

It is clear from the Figure El that the analytic approximation breaks down when we 
consider thick crystals. This is when the argument of the sincx function becomes large, the 
wave tails of the sincx cease to be insignificant - and the Gaussian approximation is not 
legitimate anymore. However for thin crystals all the analytic approximations work well. 
We will show later that this region is the most interesting when one intends to maximize 
the coupling efficiency. 

In the FigurelUwe have shown the coincidence probability as a function of the pump beam 
waist Wp. Solid lines again represent the results of the analytical approximation (fTH|l . dashed 
lines - the analogous numerical results. The plots demonstrate several cases: the bottom 
curves - for L = 1mm and w = 50fim, the middle curves - for L = 1mm and w = 150/im 
and the top, bold curves show the probability p as a function of wp with optimized values 
of L and w. 

The figures reveal a good agreement between the analytical theory and the numerical 
calculations. The discrepancies are more relevant only in the region of thick crystals or a 
very intense focusing of the beams. 

To find the optimum choice of L and w we have numerically maximized the probability 
p given by the expression (fTH|l in the region depicted in Figure |3] and performed a linear fit 
to the obtained dependencies. As a result we have obtained: 



Wopt ~ 3.32/im + 0.308wp 
Lopt ~ 167/im + 14.8wp. 

(19) 

The bolder top curves on the Figure 0] (analytical and numerical) involve an application of 
the above dependencies. 
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From the Figure |31 one can notice that the couphng process is the most effective for 
a very small pump beam mode waist. According to the numerical result depicted in this 
figure we find that the optimum waist wp^ opt ~ 27/im which, according to the expressions 
(|19|). corresponds to Wopt = 12yum and Lopt = 0.57mm. Rayleigh range for the optimally 
focused pumping beam equals 98mm. Therefore the condition of the large Rayleigh range 
in comparison with the crystal length is well preserved. 

IV. DOUBLE CRYSTAL 

So far we have considered the case of a single nonlinear crystal. However, to generate 
polarization entanglement with the use of type-I process in the scheme proposed by Kwiat 
et al. one needs two crystals rotated in respect to each other by the angle of 90° around 
the z axis. In this scheme a pump beam is linearly polarized so that the chance for a 
down-conversion in the first and in the second crystal are exactly the same. 

If the first crystal is oriented as discussed in the previous sections (the optical axis of 
the crystal lying in a plane oriented at the angle 45° in respect to the horizontal plane), 
then we have two options concerning the orientation of the second crystal. It may be either 
horizontally or vertically flipped in respect to the geometry of the first crystal. Let us notice 
that in the latter case, when the crystals are symmetric in respect to the horizontal plane, 
the amplitude ^^{uJsi ^i) given by equation is exactly the same for both crystals. Therefore 
the distinguishability of the photon pairs from the both crystals may only be a consequence 
of either spatial separation of the emission cones or propagation effects: the photons from 
the first crystal propagate through the second one as the extraordinary rays and undergo a 
transversal shift due to the walk-ofT effect. Secondly, the pump beam polarization component 
down-converting in the second crystal and propagating through the first one is delayed in 
respect to the other polarization component. Consequently, photons generated in both 
crystals arrive to the fibers at different times. However, these two dominant propagation 
effects are related to linear phase shifts easily compensable with additional linear optics. 
The pump beam delay may be compensated by a birefringent material inserted in front 
of the double crystal, and the walk-off effect can be compensated by a pair of inverted 
down-conversion crystals placed in front of the fibers. 

At this point our motivation concerning the diagonal orientation of the crystals becomes 
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clear. In the proposed configuration only the spatial separation of the photon emission cones 
can give rise to the photon distinguishability. However for thin crystals separated by a small 
distance, one can attain the perfect indistinguishability of photon pairs emerging from the 
both crystals and consequently maximize the polarization entanglement. 

V. CONCLUSIONS 

We have developed the approximate analytic theory of a type-I down-conversion to an- 
alyze the properties of the fiber coupled photon pairs. We have shown that in the regime 
of the weak focusing one can separate the coincidence spectrum, while for strong focusing 
one can optimize the coincidence efficiency. We have analytically derived the coincidence 
spectrum separability condition and the pair fiber-coupling probability and compared the 
result with the numerical calculation obtaining a good agreement. 
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APPENDIX 

Let us introduce a following coordinate system: Z axis perpendicular to the crystal and 
transverse coordinates' axes X and F, so that the optical axis lies in the plane XZ, at the 
angle a respect to the Z axis. The dispersion relations for the ordinary and extraordinary 
beams have the form, respectively llO|: 
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ion 1 \ 

k"z{kx,kY,uj) = \l ^ kj^-k^ 



1 ■ / 1 n'=(a;)2 

sm a cos a 1 

sm a H — rf^a cos^ a 

sm a + ^^7^2 cos^ a 

where n°{uj) and n*'(c(j) are the ordinary and extraordinary indices of refraction, respectively, 
given by the Sellmeier formulas [l^ . 

From the above relations one can calculate the quantities appearing in Equation (|7j): walk- 
off angle 7, group velocity combinations A/3-i- ^ and angle of propagation of the degenerate 
photons inside the crystal 9^ defined as: 



^ ^ afc|(o, 0,2^0) 



dk 



X 



gfc|(0,0,2^o) ^gA:|(0,0,a;o) 

^H±,z = 7{ ± 



du did 



O'o = -tH- (21) 



[1] P. G. Kwiat, K. Mattle, H. Weinfurter, A. Zeilinger, A. V. Sergienko, and Y. Shih Phys. Rev. 

Lett. 75, 4337 (1995). 
[2] C. H. Monken, P. H. Souto Ribeiro, and S. Pdua Phys. Rev. A 57, R2267 (1998). 
[3] C. Kurtsiefer, M. Oberparleiter and H. Weinfurter, Phys. Rev. A 64, 023802 (2001). 
[4] F. A. Bovino, P. Varisco, A. M. Colla, G. Castagnoli, G. Di Giuseppe and A. V. Sergienko, 

|quant-ph/0303126| (2003). 



[5] S. Castelletto, L P. Degiovanni, M. Ware and A. Migdall, quant-ph/Q311Q99 (2003). 



[6] K. Banaszek, A. Dragan, W. Wasilewski and C. Radzewicz, Phys. Rev. Lett, (in press). 

14 



[7] A. B. U'Ren, K. Banaszek and 1. Walmsley, Quant. Inf. Comp. 3, 480 (2003). 
[8] T. E. Keller and M. H. Rubin Ptiys. Rev. A 56, 1534 (1997). 
[9] M. H. Rubin Phys. Rev. A 54, 5349 (1996). 
[10] A. Yariv, Quantum Electronics, 3rd edition, Wiley, New York (1988). 

[11] D. N. Klyshko, Photons and Nonlinear Optics, Gordon and Breach Science Publish., New 
York (1988). 

[12] D. N. Nikogosyan, Poperties of optical and laser-related materials, Wiley, New York (1997). 
[13] P. G. Kwiat, E. Waks, A. G. White, I. Appelbaum and P. H. Eberhard, Phys. Rev. A 60, 
R773 (1999). 



15 




FIG. 1: Scheme of the setup: a pump beam Ep{r,t) focused on the nonhnear crystal (length L) 
with a waist wp and the coupled modes with waists w and a shift h at the output face of the crystal 
in respect to the pump beam. The coupled modes s and i are oriented at the angles 9s and 6i with 
respect to the z axis and the coupling lenses are preceded with interference filters. 
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FIG. 2: Coincidence spectrum \ilj{uJs,(jJi)\'^ for several mode waist diameters. For strongly focused 
beams, the spectrum becomes inseparable, and the separability can be achieved only for unfocused 
beams. 
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FIG. 3: Relative coincidence probability as a function of the crystal length L for fixed w = wp = 
100/zm. Solid line - analytical approximation, dashed line - numerical result. The Rayleigh range 
for the pumping beam is equal to 13cm. 
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FIG. 4: Relative coincidence probability as a function of the pump beam waist wp: solid line - 
analytical approximation, dashed line - numerical result for several cases. The bottom curves- for 
L = 1mm and w = 50//m, the middle curves - for L = 1mm and w = 150/im and the top, bold 
curves for L and w chosen to maximize the result for each wp. 
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